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ABSTRACT 



A recent re-analysis of Einstein data, and new i2 T observations, have re- 
vealed the presence of at least two components in the X-ray spectra of X-ray faint 
early-type galaxies: a relatively hard component [kT > 1.5 keV), and a very soft 
component [kT ~ 0.2 — 0.3 keV). In this paper we address the problem of the nature 
of the very soft component, and whether it can be due to a hot interstellar medium 
(ISM), or is most likely originated by the collective emission of very soft stellar 
sources. To this purpose, hydrodynamical evolutionary sequences for the secular 
behavior of gas flows in ellipticals have been performed, varying the type la super- 
novae rate of explosion, and the dark matter amount and distribution. The results 
are compared with the observational X-ray data: the average Einstein spectrum 
for 6 X-ray faint early-type galaxies (among which are NGC4365 and NGC4697), 
and the spectrum obtained by the ROSAT pointed observation of NGC4365. The 
very soft component could be entirely explained with a hot ISM only in galaxies 
such as NGC4697, i.e. when the depth of the potential well - on which the average 
ISM temperature strongly depends - is quite shallow; in NGC4365 a diffuse hot 
ISM would have a temperature larger than that of the very soft component, be- 
cause of the deeper potential well. So, in NGC4365 the softest contribution to the 
X-ray emission comes certainly from stellar sources. As stellar soft X-ray emitters, 
we consider late-type stellar coronae, supersoft sources such as those discovered 
by ROSAT in the Magellanic Clouds and M31, and RS CVn systems. All these 
candidates can be substantial contributors to the very soft emission, though none 
of them, taken separately, plausibly accounts entirely for its properties. 

We finally present a model for the X-ray emission of NGC4365, to reproduce 
in detail the results of the ROSAT pointed observation, including PSPC spectrum 
and radial surface brightness distribution. The present data may suggest that the 
X-ray surface brightness is more extended than the optical profile. In this case, 
a straightforward explanation in terms of stellar sources could be not satisfactory. 
The available data can be better explained with three different contributions: a 
very soft component of stellar origin, a hard component from X-ray binaries, and 
a ~0.6 keV hot ISM. The latter can explain the extended X-ray surface brightness 
profile, if the galaxy has a dark-to-luminous mass ratio of 9, with the dark matter 
very broadly distributed, and a SNIa explosion rate of ~0.6 the Tamman rate. 
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1. Introduction 



After the launch of the Einstein sateUite (Giacconi et al. 1979), it was first 
realized through X-ray measurements that elliptical galaxies may retain a large 
amount (~ 10® - 10^^ Mq) of hot (T ~ 10^ K) interstellar gas (e.g. Forman, 
Jones, and Tucker 1985; Canizares, Fabbiano, and Trinchieri 1987; see the review 
of Fabbiano 1989). This conclusion was supported by the presence of X-ray emission 
displaced from the optical body in some cases, by the relatively soft (~1 keV) X-ray 
spectra of X-ray bright ellipticals, and by the steep correlation between X-ray and 
optical luminosity (ixociB^'^"^); see Fig. 1. In contrast in spiral galaxies, where 
the X-ray emission can be attributed to a population of evolved stellar sources 
(see Fabbiano 1989, and references therein), the X-ray spectra tend to be harder, 
and LxocLb- Besides suggesting a steep relationship between Lx and Lb, the 
Lx~Lb diagram of early-type galaxies shows a substantial amount of scatter: the 
range of Lx for a given optical luminosity spans about two orders of magnitude, 
for Lb> 3 X W^^Lq (Fig. 1). 

Many theoretical models were developed to explain the Lx-Lb diagram, in- 
cluding numerical simulations for the behavior of gas flows fed by stellar mass loss 
and heated to X-ray temperatures by type la supernovae (SNIa; see the review of 
Fabbiano 1989; Loewenstein and Mathews 1987; Sarazin and White 1987, 1988; 
D'Ercole et al. 1989; David et al. 1991; Ciotti et al. 1991, hereafter CDPR). The 
most successful models are those explaining the scatter in the Lx-Lb diagram in 
terms of different dynamical phases for the hot gas flows, ranging from winds to 
inflows (CDPR). These authors, assuming that the SNIa explosion rate is declining 
with time slightly faster than the rate with which mass is lost by stars, find that 
in the beginning the energy released by SNIa's can drive the gas out of the galax- 
ies through a supersonic wind. As the collective SNIa's energy input decreases, a 
subsonic outflow takes place, and this gradually slows down until a central cooling 
catastrophe leads to the onset of an inflow. An attractive feature of this scenario 
is that any of the three phases wind/outflow/inflow can be found at the present 
epoch, depending only on the various depths and shapes of the potential wells of 
the galaxies: in X-ray faint ellipticals the gas is still in the wind phase, and the 
emission is mostly accounted for by stellar sources; in the bulk of galaxies the gas 
flows are in the outflow phase; in X-ray brightest galaxies the soft X-ray emitting 
gas dominates the emission, being in the inflow phase. In this latter case, the gas 
flow resembles a "cooling flow" (e.g. Sarazin and White 1987, 1988). 

The knowledge of the actual contribution of discrete sources to the total X- 
ray emission is an important tool for constraining the gas flow phase, especially 
in X-ray faint galaxies. Unfortunately, even the amount of the hard contribution 
coming from accreting low-mass binaries, which should be substantial, is still a 
matter of debate, due to the impossibility of estimating it directly from the available 
Einstein X-ray spectra (Forman, Jones, and Tucker 1985; Canizares, Fabbiano, and 
Trinchieri 1987). Progress in this field comes from a recent systematic investigation 
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of the X-ray spectra of early-type galaxies obtained by the Imaging Proportional 
Counter (IPC) on board of the Einstein satellite, by Kim, Fabbiano, and Trinchieri 
1992 (see Fig. 1). These authors find that on the average the X-ray emission 
temperature increases with decreasing X-ray to optical ratio (Lx/Lb), until the 
dominant contribution to the total emission comes from a hard component, similar 
to that dominating the emission in spirals, which can be fitted with a kT > 3 keV 
bremsstrahlung model. These findings are in agreement with the CDPR scenario, 
in which the emission of X-ray faint early-type galaxies is due mainly to a stellar 
component (these galaxies would be in the wind phase), while the gaseous emission 
becomes dominant in X-ray bright galaxies, which would be in the inflow phase. 

In those galaxies with the lowest values for Lx/Lb (hereafter Group 1, which 
includes 6 early-type galaxies with 29.3 < \og[Lx/LB) < 30.0, with Lx in erg s~^. 
Lb in -£©), Kim, Fabbiano, and Trinchieri 1992 find that a very soft [kT 0.2 keV) 
thermal X-ray component (hereafter VSC) is present in addition to the hard ther- 
mal component. The VSC emission has been estimated to amount approximately 
to one third to half (~ lO"^" erg s~-^) of the total X-ray luminosity in the 0.2-3.5 
keV band. With the available data, it remains possible that this VSC is present in 
all early-type galaxies. Recent ROSAT PSPC observations of two galaxies of Group 
1 (NGC4365 and NGC4382, Fabbiano et al. 1994, the companion paper) confirm 
the existence of a VSC. 

Investigating the nature of the VSC is the principal aim of this paper. We 
present hydrodynamical evolutionary sequences for the secular behavior of gas flows 
- in the frame of the CDPR models - specific to the Group 1 galaxies NGC4697 and 
NGC4365, to explore whether such flows can account for the observed VSC char- 
acteristics. These two galaxies were chosen because of their different central stellar 
velocity dispersions ; is proportional to the central depth of the galaxy poten- 
tial well, on which the gas temperature strongly depends. In our models we adopt 
the observed optical characteristics of the galaxies, and we vary the SNIa's rate of 
explosion, and the amount and distribution of dark matter. The ROSAT results 
on NGC4365 allow us to make more secure conjectures on the various components 
of the X-ray emission, considering both the X-ray spectrum and the X-ray sur- 
face brightness profile. We also consider whether the VSC could originate in the 
collective emission of discrete very soft sources, such as late-type stellar coronae, su- 
persoft accreting binaries like those detected by ROSAT in the Magellanic Clouds 
and M31 (Greiner et al. 1991; Kahabka and Pietsch 1992; Schaeidt et al. 1993), 
and RS CVn systems. 

This paper is organized as follows: in §2 we summarize and re-examine the 
observational properties of the VSC in X-ray faint early-type galaxies; in §3 we in- 
vestigate whether these properties are consistent with the expected characteristics 
of gas flows in Group 1 galaxies; in §4 we show the results of the numerical sim- 
ulations to model the observations of NGC4697 and NGC4365; in §5 we consider 
stellar sources as origin of the VSC; in §6 we present a detailed model for the X-ray 
emission of NGC4365; in §7 we summarize our conclusions. 
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2. The VSC from Einstein and i2 05^ T Observations 



For an accurate theoretical investigation, we need to have the best definition 
of the VSC characteristics which is possible within the observational boundaries. 
So, our first step was to re-analyze in detail the average spectrum obtained from 
Einstein IPC data for Group 1 galaxies by Kim, Fabbiano, and Trinchieri 1992, to 
determine the range of the VSC temperature (Tysc) and the VSC contribution to 
the total X-ray emission {Lx,vsc) allowed by the data. To this purpose, we used 
the PROS ('xray') package oilRAF, developed at SAO for the analysis of reduced 
X-ray data. Following Kim, Fabbiano, and Trinchieri 1992, we adopt an average 
value for the hydrogen column density along the line of sight Nh = 3.2 X 10^" cm~^. 
A one-temperature optically thin thermal model with solar abundance gives an 
unacceptable best-fit %^ (minimum %^ = 17.5 for 7 degrees of freedom). Acceptable 
fits are obtained when the spectrum is modeled by the superposition of a hard 
thermal component [kT > 1 keV), and a soft thermal component [kT < 0.5 keV), 
for solar abundance of the emitting gas; in Fig. 2 the best fit case is represented 
(minimum %^ = 7.9 for 5 degrees of freedom)^. The ranges of values given above 
for the temperatures correspond to the 90% confidence level, for three interesting 
parameters: the two temperatures, and the ratio of the normalizations of the two 
thermal components (Fig. 3). The temperatures are not well constrained for very 
low and high values, because the combined mirror and IPC spectral response drops 
significantly below 0.5 keV and beyond 4 keV, and the IPC has very poor spectral 
resolution. 

The uncertainties in Group 1 data are such that the normalizations of the two 
thermal components are not well determined. As an example, if we keep 1 < kT < 5 
for the hard component, and 0.1 < kTysc < 0.5 (see Fig. 3), we have that 
0.15 < Lx,vsc / Lx < 0.50 at the 90% confidence level. Assuming that this range 
is representative of all Group 1 galaxies, we derive the range of Lx,vsc values given 
in Table 1 for NGC4697. 

Recent i2 05^ T observations of two Group 1 galaxies (NGC4382 and NGC4365) 
confirm the presence of a VSC with better spectral data. We do not describe these 
observations here, because they are discussed in detail in the companion paper 
(Fabbiano et al. 1994). Again a one-temperature model is excluded, except for the 
case of a temperature between 0.6 and 1.1 keV, and abundance values close to zero. 
Another more realistic possibility (Fabbiano et al. 1994) is the combination of a 
hard [kT > 1.5 keV) and a soft thermal component, with kTysc = 0.2 keV at 
the best fit. In Fig. 4a we show the ranges of variation at 68%, 90% and 99% 
confidence level for kTysc in the NGC4365 case (line of sight absorption has been 
assumed, see Table 1). A comparison with the similar figure obtained from the 

^ With these data we cannot model temperatures and metal abundance of the 
ISM meaningfully. However the ROSAT observations of two Group 1 galaxies 
(Fabbiano et al. 1994) suggest that the two temperature model is the more likely 
representation of the data. 
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Einstein IPC data (Fig. 3), shows how the contours in Fig. 4a are still open 
towards high temperatures, due to the PSPC spectral response, but now the VSC 
temperature is well constrained. The soft and the hard contributions to the total 
X-ray luminosity in the ROSAT band are comparable in the best fit case, but can 
vary up to a factor of two in each direction (Fabbiano et al. 1994). The spectral 
counts distribution of NGC4365 together with the two best fit components and their 
sum are shown in Fig. 4b. We report in Table 1 the main results of the analysis 
concerning NGC4365, for which we perform numerical simulations of the hot gas 
behavior (we do not simulate the SO galaxy NGC4382, because our spherically 
symmetric models cannot handle the presence of a disk). 

3. Can a hot ISM be the origin of the VSC? 

In this section we examine whether the observed VSC characteristics are in 
principle compatible with a hot ISM origin. As an outcome of this investigation, 
we will also find what input parameters are needed for the numerical simulations 
of the gas flow behavior. We first introduce our galaxy models. 

3.1 The galaxy models 

The galaxy structure adopted for the hydrodynamical simulations of the hot 
gas evolution is that of CDPR: it is a superposition of a King 1972 model for 
the luminous matter distribution, and of a quasi-isothermal halo for that of the 
dark matter; both distributions are truncated at the same tidal radius rt. We 
adopt rt/rc*=180, where r^* is the stellar core radius, so that the projection of 
the luminous matter distribution gives a good fit to the observed optical surface 
brightness profiles of ellipticals (Kormendy 1982). The observed quantities are 

ct*; and rc*; from these the central stellar density /jq* is derived through the 
virial condition (e.g. Binney and Tremaine 1987); in this way the luminous matter 
distribution is completely defined, because the stellar mass-to-light ratio is assumed 
to be constant with radius. The dark matter distribution can be determined by 
choosing the ratios of the dark/luminous masses (i2 = Mh/M^), the dark/luminous 
core radii (/3 = Vch/fc*), and the dark/luminous central densities (7 = poh/po*), 
only two among these three parameters being independent. The problem of the 
physical consistency of the proposed mass distribution - whether it represents a 
self-consistent solution of the Poisson and Vlasov equations (see, e.g., Binney and 
Tremaine 1987) with two components - was addressed by Ciotti and Pellegrini 1992, 
in the hypothesis of stellar orbits with various anisotropies. In this paper we just 
want to derive the physical properties of the gas flows, and possible constraints on 
the amount and distribution of the dark matter, starting from reasonable shapes 
for the potential wells. 

The time-evolving input ingredients of the numerical simulations are the rate 
of stellar mass losses, and the SNIa's heating; these too are calculated as in CD- 
PR. In particular, the stellar mass loss rate is parameterized as M^[t) ~ —1.5 X 
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10 ^^L^it/ 15 Gyvs) -^'^ Mqjv ^] the SNIa's heating is parametrized as LsN{t) — 
7.1xlO^°'i?SNiB(Vl5Gyrs)-^-s erg s'^ with < i^sN < 1- The value t?sn = 1 gives 
the present rate of SNIa explosion in ellipticals estimated by Tammann 1982; this 
is O.SS/i^SNU, where h = ITo/lOO, and 1 SNU=1 event/100 yrs/10^° Lq (we adopt 
Ho = 50 km/sec/Mpc throughout this paper). Van den Bergh and Tammann 1991 
give an estimate which translates into i^sN = 1-1; close to that of Tammann 1982; 
Cappellaro et al. 1993 give a value 4 times lower (i^sN ~ 0.3), which illustrates the 
present range of uncertainty in the SNIa rate of explosion in ellipticals. Lsn is 
assumed here to be decreasing with time just faster than the stellar mass loss rate 
(Renzini 1989). This assumption produces a time decreasing specific heating of the 
gas, and so the wind/outflow/inflow secular evolution of the gas dynamical state, 
described in §1. 

The effects introduced by the possible presence of population gradients can 
be neglected without affecting our conclusions. Age and/or metallicity gradients - 
such as those that can be estimated from observations - introduce a variation of 
M^{f) within ~ 10% (Renzini and Buzzoni 1986), while many uncertainties already 
affect LsNif)- 

In summary, our parameter space for each galaxy model consists of (i2, /3, i^sn); 
the first two parameters describing the dark matter amount and concentration, the 
third describing the present rate of SNIa explosions. 

3-2 Lx,gas and Lx,vsc 

Here we address the question of which dynamical phase of the gas flow can 
reproduce the observed Lx,vsc- Global inflows are far too luminous (-£x,£/as~10'^^ 
erg s~-^) with respect to Lx,vsc, while during winds the X-ray luminosity of the 
gas is very low (< 10^^ erg s~-^); during the intermediate outflow phase Lx,gas 
can vary by orders of magnitude, going from wind values to inflow values. The 
duration of this transient phase is very sensitive to the input parameters, and it 
can last up to ~10 Gyrs; this means that any of the Lx values, from those typical 
of winds to those typical of inflows, can be shown by gas flows in outflows at the 
present epoch, depending on the input parameters. Another dynamical state in 
which Lx,gas ~ Lx^vsc is the "partial wind" (hereafter PW). In this case radiative 
losses suppress an outflow in the center of the galaxy, causing a small inflow region, 
but a wind can be sustained in the external regions, where the gas density is much 
lower and the gas is also less tightly bound (MacDonald and Bailey 1981; §4.1.1). 
Clearly, the smaller the central inflow region, and the higher the wind velocity, the 
lower is the gas mass inside the galaxy, and so the X-ray luminosity. 

Which gas flow phase is expected in Group 1 galaxies? They show a large 
range of values for the central stellar velocity dispersion cr^ (cr^ =186-284 kms"-^; 
Whitmore, McElroy, and Tonry 1985), reaching particularly low values with re- 
spect to the average Lb-c^ relation for early-type galaxies (Faber &: Jackson 1976, 
Davies et al. 1983). This is important because is one of the key parameters 
regulating the gas flow phase; together with the amount of dark matter (-R), and 
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its concentration (/3), it measures the binding energy of the gas. Following CD- 
PR, we can summarize the gas flow energetic balance as a function of cr^, and the 
introduced parameters [R, /3, i^sn); by using the parameter %. This is deflned as 
the ratio of the power required to steadily extract the gas shed by stars from the 
potential well of the galaxy, to the energy made available by SNIa's explosions: 
X oc crl[A + /3~*' "^'^-R)/'i?sN; at any epoch [A is a constant). When % >1, inflows 
are expected, while when % < 1, we expect winds; if % 1 we have outflows. In 
Table 2 we show the values assumed by % at the present epoch, for i^sN = 1 and 
for the values of the two galaxies in Table 1, NGC4697 and NGC4365; the dark 
matter amount (R) and concentration (/3) have been varied. We can see how winds 
are expected in both galaxies, if i^sN = 1; even for large R values. So, we have to 
lower i^sN to obtain outflows, especially for NGC4697, which has the smaller a^. 

The parameter % cannot be used to make predictions on the PW phase: the 
gas could be in PW with % either larger or smaller than 1, because this parameter 
describes the global energetic balance of the gas, while in PW the gas behaves 
differently in the center and in the outskirt of the galaxy (these regions become 
"decoupled"). We expect PW solutions for small amounts of dark matter - since 
we assume it is distributed mainly at large radii - and this is conflrmed by the 
numerical simulations (see also §4.1.1 and 4.2.1). 

In conclusion, we have to explore through numerical simulations of the gas flow 
evolution which region of the parameter space [R, /3, i^sn) is populated by gas flow 
phases with Lx,gas ~ -^x,ysc! at the present epoch. We can just anticipate we 
need i^sN < 1- 

3.3 Tgas and Tvsc 

We explore here whether the VSC temperature is compatible with that of the 
hot ISM in ellipticals. As it was the case for the X-ray luminosity, the average gas 
flow temperature depends strongly on a^. For instance, during the inflow phase 
predictions on the temperature proflle of the gas can be made where the following 
simplifying assumptions hold: 1) the gas is in quasi-hydrostatic equilibrium in 
the galaxy potential well (i.e. tgoundif) « tcoolif), with t 

sound the sound speed 

crossing time, and tcool the gas cooling time); 2) the flow is dominated by gravity 
rather than by the supernovae energy input (the gas is in inflow); 3) the cooling 
is just a small fraction of the global energetic balance of the gas flow. Under 
these hypotheses, one derives for the gas Tin{r) oc |0(t")|, (j) being the total galaxy 
gravitational potential. Just outside the central region of the galaxy, as soon as the 
infall velocity becomes small and the above hypotheses apply, one has 
3.6fimpalF[R, (3), where fimp is the average gas particle mass, and F[R,(3) is a 
form factor of the order of unity; F{R,l3) = 0.46 ^ 1.6, as i2 = ^ 10, and (3 = 2. 
So kTin{5rc*)^0.3 keV if cr^^lSO km . 

In conclusion, we can expect the gas flows to be characterized by temperatures 
of the order of Tvsc only for low values. However the temperature estimates 
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given above refer to the central regions of the galaxy, while the gas radial temper- 
ature distribution decreases outward, unless a sizable external pressure is effective, 
or a very large amount of dark matter broadly distributed (see, e.g., Sarazin &: 
White 1987). So, a detailed calculation of the X-ray spectrum of the gas in the 
sensitivity bands of the satellites is required. The extent to which the dark matter 
presence and distribution, the SNIa's rate, and the cooling of the gas determine the 
past flow evolution - and then its present temperature and density profiles - can 
be derived only through numerical simulations. 



4. Numerical simulations for an ISM origin of the VSC 

In this section we explore the gas flow evolution in the parameter space [R, /3, 
"i^sn); to find whether there are regions in which gas flows have the observed VSC 
characteristics at the present epoch. To this purpose we construct hydrodynamical 
evolutionary sequences specific to the galaxies in Table 1, NGC4697 and NGC4365; 
these have been chosen because, in Group 1, they have respectively the lowest and 
one of the highest values for cr^, on which Lx,gas and Tgas mostly depend (§3). To 
tailor the model structures to these galaxies, we adopt also the observed values of 
the blue luminosities (Table 1) and of the core radii (King 1978). Initially only the 
two fundamental parameters R (dark-to-light mass ratio) and i^sN (SNIa's rate) 
vary, while the dark-to-luminous core radii ratio /3 = 2; we explore later the effect 
of a variation of /3. The mass profiles of such galaxy models are shown in Fig. 5, 
for /3 = 2 and /3 = 10. 

The time-dependent equations of hydrodynamics with source terms, and the 
numerical code to solve them, are described in CDPR. The cooling curve is an 
analytical fit to the results given by the Raymond code for the thermal emission 
of a hot, optically thin plasma, at the coUisional ionization equilibrium, and with 
solar abundance (Raymond and Smith 1977). The gas flow evolution is followed for 
15 Gyrs. 

The characteristics of the flows (X-ray luminosities in the Einstein and ROSAT 
sensitivity bands, spectral energy distributions, and X-ray surface brightness pro- 
files) are compared directly to the data, after they are transformed into predicted 
"observed" quantities, as follows. Using the Raymond code, the spectral energy 
distribution in the bands of sensitivity of the satellites is computed in each one 
of the 100 grid points used for the hydrodynamical simulation. After projection, 
the energy distribution of the photon flux is corrected for the interstellar absorp- 
tion along the line of sight, using the photo-absorption cross sections of Morrison 
and McCammon 1983. Finally the model spectral distribution is convolved with 
the Einstein IPC (Harnden et al. 1984) or ROSAT PSPC spectral response and 
effective area, to obtain the predicted counts distribution. The model X-ray sur- 
face brightness profile is convolved with the energy dependent PSPC point spread 
function (Hasinger et al. 1992). Such predicted models are compared directly to 
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the observational data, and the results of the comparison are described below, in 
§4.1.2, 4.2.2, and §6. 

4.1 NGC4697 
4.1.1 The X-ray Luminosity 

In Fig. 6 we show how the various gas flow phases, at the end of the simulation, 
populate the [R, i^sn) plane; the regions where Lx,gas{^^ Gyrs)~ Lx,vsc are also 
shown. Due to the low value, a large region of the plane is occupied by global 
winds; to obtain %(15Gyrs)sal, that is to expect an outflow at the present epoch, 
we need a dark-to-light mass ratio i2 40, if i^sN = 1- 

For R > Q the gas flow evolution is the "classical" one: wind/outflow/inflow; as 
anticipated in §3.2, Lx,gas{^^ Gyrs)sa Lx^vsc during the outflow. The i^sN value 
for this is very critical, because the flow is very sensitive to the global energetic 
balance. As an example, for [R=9, 'i?sN=0.38) the gas is in outflow, and its Lx can 
vary by a factor of 8 with i^sN varying by only 1%. 

If R^Q the flow doesn't follow the evolutionary history wind/outflow/inflow. If 
i^SN^O.S, i.e. excluding the cases when a global wind lasts up to the present epoch, 
a cooling catastrophe takes place in the flrst few Gyrs (Fig. 7) in the central regions 
of the galaxy, while the gas escapes the external regions at a high velocity up to 
the present, and we have a PW. This is due to the SNIa's heating being too low 
to prevent an inflow in the gas dense and tightly bound central regions, while it is 
high enough to drive out the gas in the external regions. The tendency is however 
to evolve into a subsonic outflow, and then eventually into an inflow, even for the 
outer regions, as the speciflc heating secularly declines (§3.1). 

The flow behavior across the [R, i^sn) plane can be explained looking at how 
the structure of the models changes across it. In the central regions of the galaxy 
the amount of luminous matter is always larger than that of the dark matter (see 
Fig. 5); increasing R, we bring to inner radii the dominance of the dark matter 
over the luminous one, but even for R=9, we have Mh{< r)/M^[< r) > 1 only for 
r > lire*. So, the central regions are always dominated (from a gravitational point 
of view) by the stellar structure, which is flxed across the plane. At low values of 
R, to avoid winds we decrease i^sN; until it reaches a critical value i^p^y which is 
too low to prevent an inflow in the central regions, and we have a partial wind. At 
high values of R, T^p^ cannot prevent a global inflow, because also the external 
regions become very tightly bound. 

The above results are sensitive to variations of the parameter /3 especially 
during the outflow phase: a larger value for /3 - decreasing the dark matter con- 
centration: Mh{< r)/M^{< r) > 1 for r > 21rc*, if = 9 and /3 = 10 (Fig. 5) 
- has the effect of retarding the onset of the global inflow. As a consequence, the 
values of i^sN corresponding to outflows should be slightly decreased in Fig. 6. In 
PW cases - when R is low, and the luminous matter concentration is already high 
enough to cause a central inflow - a /3 > 2 has little influence. Values of /3 smaller 
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than ~ 2 are not plausible, since dark matter halos around galaxies (if present) 
seem to be diffuse (Mathews 1988; Bertin, Saglia, and Stiavelli 1992). 

4.1.2 The X-ray Spectrum 

The predicted pulseheight distributions of a few representative cases of gas 
flows having Lx,gas{^^ Gyrs)sa Lx,vsc are shown in Fig. 8. Here also the pulse- 
height distribution corresponding to the best fit temperature of the VSC from 
Einstein IPC data is shown, together with that corresponding to the 90% upper 
limit on Tysc- The predicted distributions are closer to the best fit distribution for 
those models with low values of i^sN and R, i.e. PWs, when the gas temperature 
is lower. For i2>5 the peak in the counts distribution of the gas flow models takes 
place at energies higher than in the best fit case, and it is clearly displaced for R = 9, 
"i^SN = 0.38. However, all the predicted distributions lie within that corresponding 
to kT = 0.5 keV, the 90% confidence upper limit on the VSC temperature, because 
the average X-ray emission temperature of gas flows is kT = 0.2 — 0.3 keV. In Fig. 9 
we show the comparison of the average Group 1 spectrum with that of the [R = 3, 
"i^SN = 0.2) model, to which a hard thermal component has been added. 

In conclusion, for a galaxy with the optical characteristics of NGC4697, an ISM 
origin of the VSC is in principle possible; the strongest constraints in the param- 
eter space come from the required Lx,gas value, not from that of the gas emission 
temperature. In fact the statistical uncertainties in deriving the temperature of the 
VSC from Einstein IPC data are comparable to the spread of the model tempera- 
tures. So, these large uncertainties don't allow any discrimination between different 
values of R and i^sN- 

4.2 NGC4365 

This galaxy, belonging to the Virgo cluster, has been observed also by the 
ROSAT PSPC. In the modelling we try to reproduce the results of these more 
accurate observations. 

4.2.1 The X-ray Luminosity 

In Fig. 10 we show the (i2, i^sn) plane for NGC4365. The main features of 
the gas flow behavior are the same as in the case of NGC4697; here we point out 
the main differences, and some characteristics of the gas flows which have been 
investigated in more detail, since more accurate observations are available. 

Due to the larger value of cr^, a larger region of the (i2, i^sn) plane is occupied 
by global inflows, and even with i2 ~ the flows can be luminous enough (or 
even too luminous); the t^p^y values are higher than for NGC4697; the models 
with ix,£/as(15 Gyrs)~ Lx,vsc are in PW up to i2 5. For 3^-R < 5 the global 
wind phase is actually followed by a global outflow phase, but the emission of the 
latter is never sufficiently high: a central cooling catastrophe develops, originating 
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a PW, before the outflow phase can reach Lx,vsc (Fig. 11). Global outflows 
can have -£x,£/as~6 X 10^^ erg only for i2>5. The i^sN value to have Lx,gas{^5 
Gyrs)~ Lx,vsc is specially critical ii R = 5 — 6, the lowest R values for which 
the outflow can be luminous enough: enough gas mass has accumulated within the 
galaxy only just before the central cooling catastrophe occurs (see the R = 5 case 
in Fig. 11). 

As for NGC4697, changing the value of /3 has some effect only during outflows; 
for instance, ii R = 9 and we increase /3 from 2 to 10, i^sN goes from 0.68 to 0.59; 
if -R = 6, increasing /3 from 2 to 20 lowers i^sN from 0.55 to 0.43. 

4-2.2 The X-ray Spectrum 

In Table 3 we show the average emission temperatures for a few representative 
models having Lx,gas{^^ Gyrs) ~ Lx,vsc in the i^O^^T band; the emission region 
considered is the same used by Fabbiano et al. 1994 to derive the observed spectral 
counts distribution (i.e. a circle of 250" radius). Going from PWs to outflows, kT 
varies from 0.4 to 0.7 keV. 

In Fig. 12a we plot the observed counts distribution, together with the super- 
position of the spectrum of the gas flow - corresponding to a galaxy model with 
R = Z and i^sN = 0.35 - and the spectrum of a bremsstrahlung component with 
kT ~ 4 keV, both convolved with the ROSAT response. The interstellar absorption 
is flxed at the line of sight value given in Table 1; the temperature and normal- 
ization of the hard component have been chosen so as to minimize the %^ of the 
flt with the observed counts distribution. From Fig. 12a we see how difficult it is 
to reproduce the observed peak at energies around 0.2 keV; this is true for every 
galaxy model. The case in Fig. 12a can be brought into much better agreement 
with the observational data if a column density Nh = 3 X 10"^^ cm~^ is assumed 
(Fig. 12b). 

In conclusion, the gas spectrum is not as soft as that of the VSC, and cannot 
explain the observed counts distribution at low energies, under reasonable condi- 
tions, even when no dark matter is present, and with a very low SNIa rate (see the 
R = Q and i^sN = 0.18 case in Table 3), because of the value of a^. 

4-3 Conclusions 

The numerical simulations we have performed, tailored precisely on a couple 
of galaxies in which the VSC has been measured, show that hot gas flows can 
easily be characterized by a luminosity Lx,gas ~ Lx,vsc, for acceptable values 
(within the present uncertainties) of the input parameters describing the amount 
and distribution of dark matter, and the SNIa explosion rate. The analysis of the 
X-ray spectra of gas flows reveals instead that these flows can be responsible for 
the VSC only in those galaxies having a low value of cr^ (<200 km s~-^), such as 
NGC4697. Among these galaxies, the requirement Lx,gas ~ Lx,vsc implies a low 
value of T?sN (^0.4 if R<10). 
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ROSAT observations, confirming the existence of a VSC, indicate that its 
temperature doesn't depend on the value of the host galaxy: at 99% confidence 
level kTvsc < 0.3 keV in NGC4365, whose value prevents the average gas 
emission temperature from falling below 0.4 keV. This fact leads us to conclude 
that by itself a hot ISM cannot generally explain the VSC. 

5. Stellar Sources as Origin of the VSC 

In the following we explore the alternative possibility that discrete sources 
have the required spectral characteristics and luminosity to account for the very 
soft emission in early-type galaxies. ROSAT has observed supersoft sources in 
M31 and the Magellanic Clouds, and the presence of a soft component in RS CVn 
systems has been well established, while late type stellar coronae have been known 
to be sources of very soft X-ray emission since the Einstein observations. 

5.1 Supersoft Sources 

i2 T observations, during the All-Sky Survey and deep pointings, have re- 
vealed strong supersoft X-ray emission in several stellar sources; in a fraction of 
them such emission had been revealed previously by Einstein observations. In 
the Large Magellanic Cloud (LMC) the supersoft sources are CAL 83 and CAL 
87, optically identified with close binary systems, and RXJ 0527.8-6954 (Greiner 
et al. 1991), and the variable RXJ0513.9-6951 (Schaeidt et al. 1993). In the Small 
Magellanic Cloud (SMC) supersoft sources are PN67 (at the location of a plan- 
etary nebula), 1E0035.4-7230, and SMC3, a symbiotic star (Wang and Wu 1992; 
Kahabka and Pietsch 1992). In our Galaxy ROSAT detected GQ Mus, a recent 
nova (Ogelman et al. 1993). These are among the softest X-ray sources observed, 
and it has been suggested that they form a separate new class of X-ray object- 
s, not previously encountered in our Galaxy because of the interstellar extinction 
(Schaeidt et al. 1993). Their similar spectra - the bulk of the X-ray emission is 
below 0.5 keV - can be described by black body radiation with kTbi, = 20 — 50 
eV and X-ray luminosity from a few X 10^^ erg up to ~ 10^^ erg . Finally, 
15 supersoft sources have also been discovered in M31 (G. Hasinger 1993, private 
communication) . 

The origin of the soft spectra and the nature of the compact object in such 
supersoft sources are still unclear: they could be neutron stars accreting matter 
at or above the Eddington rate (Greiner et al. 1991), or they might represent a 
long sought class of white dwarfs (WDs) burning accreted matter (van den Heuvel 
et al. 1992). This second interpretation seems more likely, thanks also to the nature 
of the accreting object optically identified in a few cases (see for example the case of 
SMC3 or that of the nova GQ Mus). In the van den Heuvel et al. model for CAL83 
and CAL87 the X-ray emission is produced by steady nuclear burning on the surface 
of a 0.7 1.2 Mq WD; the accreted matter is not ejected in thermonuclear nova 
flashes, thanks to the high accretion rate [(l-4)xlO~^M0yr~-^]. 
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5.1.1 Can Supersoft Sources Be the VSC Origin? 



We examine here whether the presence of supersoft sources can explain the 
X-ray spectra of those early-type galaxies where the VSC was discovered. Fitting 
the average Group 1 X-ray spectrum with a black body component of kTbi, = 50 
eV, instead of a 0.2 keV thermal component, the goodness of the fit remains the 
same; kTbi, = 150 eV at the best fit. The contribution of the black body component 
is approximately half of the total flux in the (0.2-4) keV band, and so, for example, 
500 supersoft sources are needed to explain the soft emission of NGC4697. 

Fitting the X-ray spectrum of NGC4365 with a black body component of 
kTbb ~ 50 eV, in place of a very soft thermal component, the fit gets worse: if 
we fix Nh at the line of sight value, the minimum %^ increases from 24.9 (with 18 
degrees of freedom) to 37.8 (with 19 degrees of freedom). Allowing kTbi, to vary, at 
the best fit kTbi, ~ 110 eV, and still the fit is slightly worse than with a 0.2 keV 
thermal component (minimum %^ = 26.7 with 18 degrees of freedom). In Fig. 13 
(which can be compared to Fig. 4a) the situation is summarized. If kTbi, ~ 110 eV, 
Lx,bb ~ 3 X lO"^" erg s~^, and so we need 300 supersoft sources to account for 
the soft emission of NGC4365. 

In conclusion, the supersoft sources can plausibly account only for a fraction 
of the very soft emission of early-type galaxies, because they are slightly too soft. 
Also, something more needs to be discovered concerning their nature (stellar com- 
ponents, origin and duration of the soft phase, variability) to investigate whether 
their collective emission can be significant in early-type galaxies. For example, the 
interpretation model of van den Heuvel et al. 1992 requires as donor stars respec- 
tively a normal F-G star with mass ~ 1.4 — 1.5 Mq for CAL 87, and either an 
evolved main sequence star of ~ 1.9 — 2.Mq, or a post main sequence subgiant of 
lower mass (M>1.5M0) for CAL 83. Such companions are unlikely to be present 
in the old population of an elliptical galaxy. In other cases the estimated compan- 
ion mass is much lower [see, e.g., the case of GQ Mus, a classical nova, where the 
mass of the secondary is < 0.2 Mq, Ogelman et al. 1993]. Recently, Rappaport 
et al. (1994) estimate that in our Galaxy today there should be more than 1000 
systems with the characteristics required by the van den Heuvel et al. model. 

Finally, it is interesting to note that were the supersoft sources WDs burning 
accreted matter, they would represent a long sought class of sources, which could 
also be the progenitors of type la supernovae (Munari and Renzini 1992, Kenyon 
et al. 1993). In this case the accreting WD would explode after having grown to 
the Chandrasekhar limit, or even before (helium detonation SNIa's). We estimate 
here what SNIa rate of explosion {Rexp) is expected if Lx,vsc comes entirely from 
WDs accreting matter and then exploding as SNIa's: Rexp ~ Lx,vsc /(>-Macc 

s~^, where a is the fraction of the bolometric luminosity produced by the hydrogen 
burning, collected by the ROSAT satellite; Mace is the mass which has to be 
accreted before explosion, in solar masses; and the energy yield per unit mass by 
hydrogen burning, for 0.7 hydrogen mass fraction, is ~ 4.5 X 10^^ erg g~^. On 
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average, the progenitor WD has to accrete ~ 0.2 Mq before exploding (Munari 
and Renzini 1992), or even less if it undergoes helium detonation (Woosley and 
Weaver 1993); so R^xp ^ 1-5 x lO'^V" s"^ in NGC4365. Since the observed rate 
is RsN ~ 4.5 X 10~"^*''i?sN s~"^ for this galaxy (from §3.1), to reproduce it we need 
a<0.03/i?sN- Note that a = 0.01,0.24,0.82 for kTu = 10,20,50 eV, and so most 
vsc should come from objects with quite low temperatures. This requires to 
make another check, i.e. whether the UV flux that would be produced is consistent 
with that observed. Greggio and Renzini (1990) derive for ellipticals Ljjv / Lbol < 
0.02; imposing this constraint on the UV emission of supersoft sources, we derive 
a > 2 X 10~^ in NGC4365. So, the possible link between Lx,vsc, supersoft sources, 
and SNIa's is not ruled out by current observational constraints. 

5.2 Late Type Stars 

In comparison to other types of stellar X-ray sources which are powered by 
accretion, normal stars are very faint sources of X-ray emission (-tx~ 10^^ — 10^" 
erg s~^). Their collective emission can still be substantial, because dwarfs of all 
spectral types from F to M are sources of X-ray emission (e.g. Schmitt et al. 1990); 
particularly, it is interesting to consider the X-ray properties of K and M dwarfs. 
These stars are the most numerous in a single burst old stellar population like that 
of elliptical galaxies. In such stars the X-ray emission is thermal, and originates 
in the stellar coronae. From a systematic investigation of Einstein IPC data, 
Schmitt et al. 1990 find that a two-temperature description is always required for 
the coronae of K and M dwarfs when high signal-to-noise spectra are available. The 
values of these temperatures cluster around kT ~ 1.4 keV and kT ~ 0.2 keV. 

We can estimate the collective emission due to the lower temperature compo- 
nent of K and M dwarfs in ellipticals by multiplying their expected number {Nk and 
Nm respectively) by their average X-ray luminosity (< Lx >K and < Lx >m)- 
Assuming that K and M dwarfs were generated with the Salpeter initial mass func- 
tion, their expected number in a single burst stellar population is: 

/.0.85 

Nk= I A-M-^-^^dM (1) 

Jo.6 

and 

/.0.6 

Nm= A-M-^-^^dM (2) 
t/o.os 

where A = 2.9Lb (CDPR), M = O.85M0 is the turn-off mass after 15 Gyrs, M^^/ = 
O.O8M0. So we derive Nm = 60.7-tB and Nk = I.QILb, with Lb in Lq] in the 
following we consider just the M dwarf contribution, since < Lx >m and < Lx >k 
are similar. Adopting < Lx >M= 2 X 10^^ erg s~^ (Rosner et al. 1981), we obtain 
Lx,M = Nm < L >M= 1-2 x IO^^Zb- This means Lx,m ~ 2 x 10"^^ erg s~^ for 
NG'C4697, and Lx,m ~ 8 x 10^° erg s'^ for NGC4365: ' even if just a fraction of 
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Lx,M is due to the 0.2 keV component, M dwarfs seem to be able to explain the 
VSC in early-type galaxies. In general, log[Lx,M / Lb) ~ 30, while for Group 1 
galaxies log{Lx,vsc / Lb) = 29.0 - 29.7. 

Two comments are however in order. The first is that in late-type stars the 
emission is strongly dependent on the rotation rate (Pallavicini et al. 1981), be- 
cause coronal heating is strongly related to the presence of a magnetic dynamo 
(Rosner 1991). Probably an "age effect", reducing the angular momentum, also 
reduces the X-ray luminosity, and this is particularly crucial for population II M 
dwarfs. Barbera et al. 1993 have systematically investigated the presence of this 
effect: using all the available Einstein IPC data, they have identified a subsam- 
ple representative of the population of M dwarfs in the solar neighborhood; using 
a kinematic criterion to determine the age, they have confirmed the decrease in 
X-ray luminosity going from young-disk to old-disk population. For old M dwarfs 
they find < Lx >M= (0.8 — 1) x 10^^ erg (going from late to early M). We 
derive then log(i:'^'^^/i:B) ~ 28.7; for NGC4365 L^f-^^ ~ 3.5 x 10^^ erg s'^ while 
Lx,vsc ~ 2.7 X lO"^" erg . So the age effect probably rules out the M dwarfs as 
the main origin of the VSC in ellipticals; the contribution of these stars could even 
be much lower, when we consider that the ROSAT pointed observation of an old M 
dwarf such as the Barnard star shows that its Lx = 3 X 10^^ erg (T. Fleming 
1993, private communication). 

The second consideration is that the Salpeter initial mass function probably 
doesn't hold down to stellar masses as low as 0.08 Mq. Tinney et al. 1992 find that 
in the solar neighborhood the initial mass function clearly flattens for stellar masses 
< O.2M0. If this is the case also in early-type galaxies, we have overestimated Nm- 
Unfortunately we don't really know what the IMF is in ellipticals; we can just note 
that the adopted Salpeter IMF would imply a mass-to-light ratio of M^/Lb ~ 13, 
within the range of values recently estimated by Guzman, Lucey, and Bower 1993 
(M^/Lb = 8 ^ 25, assuming Ly = 1.Q7Lb as Greggio and Renzini 1990). If we 
want to explain the VSC with old M dwarfs, we have to increase Nm by at least 
of a factor of 4 (8 in NGC4365). This implies an IMF steeper than the Salpeter 
function, with slope > 2.6 in eq. (2), which in turn implies M^/Lb > 19. 

In conclusion, late type stellar coronae do have the required spectral charac- 
teristics to explain the VSC in ellipticals, but due to an age effect we expect their 
collective emission to be too low to account for it all, unless the IMF for M dwarfs 
is much steeper than the Salpeter function, and hence is at variance with the IMF 
in the solar neighborhood. 

5.3 RS CVn systems 

RS Canum Venaticorum (RS CVn) systems are chromospherically active ob- 
jects, typically consisting of a G or K giant or subgiant, with a late type main 
sequence or subgiant companion (Linsky 1984). As Einstein observations have 
shown, RS CVn's are the most X-ray luminous late type stars {Lx ~ 10^^ — 10^"^'^ 
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erg s~^), and their quiescent coronae can be well modeled by a thermal plasma with 
two temperatures (see e.g. Schmitt et al. 1990). Recently, Dempsey al. 1993 pre- 
sented the results of their analysis of X-ray spectra of 44 RS CVn systems obtained 
during the ROSAT All-Sky Survey with the PSPC. They find that a bimodal 
distribution of temperatures centered near 0.2 and 1.4 keV fit the data best. 

We examine here the soft X-ray emission expected from RS CVn systems in 
early type galaxies. Though most of Dempsey et al. 's systems belong to population 
I, we take 4 X 10^" erg as an average value for the soft X-ray luminosity of a 
system (Dempsey et al. 's Fig. 3). The spatial density of RS CVn's has been 
estimated by Ottmann and Schmitt 1990 for the Galaxy; here we estimate it for 
an old stellar population as follows. The number of RS CVn is a fraction of the 
number of G or K giants and subgiants (respectively Ngiant and Nsubg)] Ngiant 
and Nsubg can be calculated from the knowledge of the time spent by stars in 
such evolutionary phases (see e.g. Renzini 1989). Assuming that low mass metal 
rich stars spend approximately 1 Gyr on the subgiant branch, and 1 Gyr on the red 
giant branch, we have Nsubg -\- Ngiant ~ O.I-Lb- Taking half this number for those in 
binary systems, and another factor (a) for those which become RS CVn, we derive 
the collective contribution Lx,RSCVn = 4 X lQ^^{Nsubg + Ngiant) ~ 2 X lO^^aie; 
or log {Lx,Rscvn/ Lb) = 28.6 if a = 0.2. So we conclude that RS CVn systems 
represent another class of likely contributors to the VSC in early-type galaxies, 
especially thanks to their spectral characteristics, but they are certainly not the 
unique explanation of the phenomenon. 

6. A Model for the X-ray Emission of NGC4365 

We propose here a model for the X-ray emission of NGC4365, to explain the 
results obtained from the ROSAT pointed observation of this galaxy (Fabbiano 
et al. 1994). The observed X-ray spectrum and surface brightness profile (Sx(-R)) 
allow us to make some conjectures about the dynamical phase of the gas flow, and 
the different components of the X-ray emission. From the resulting model we also 
obtain information on the dark matter content of this galaxy, its distribution, and 
the SNIa rate. 

6.1 The X-ray surface brightness profile 

Tix{R) obtained by ROSAT is shown in Fig. 14. Assuming that the hard 
component present in the X-ray spectrum (§2) is due to low-mass accreting binaries, 
its spatial distribution probably follows that of the stars; the same holds for the 
VSC, if it comes from stellar sources. In Fig. 14 we also plot the superposition 
of two brightness profiles, having each the spectral energy distribution and the 
normalization of the two components derived from the best fit to the observed 
spectrum; their spatial distribution is the same as that of the optically luminous 
matter (King 1978), and has been convolved with the PSPC response. Taken at face 
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value, the observed Sx(-R) appears less peaked than that of this model. However, 
the discrepancy is at the center, where eventual uncertainties on the PSF are more 
crucial. The way to really solve the point is to observe the galaxy with the ROSAT 
HRI, which we will propose to do. For the purpose of this paper, we will take 
this discrepancy at face value, and look for possible explanations. This exercise 
demonstrates the complexity of solutions one may have. 

The first way out comes from assuming that some cold material is present at 
the center of the galaxy to absorb the X-ray radiation (see, e.g.. White et al. 1991). 
In this case, the 'observed' Sx(-R) could look like that in Fig. 15, where the 
absorbing material is confined within a spatial radius of 4 kpc, and produces a 
column density of ~ 5 X 10^^ cm~^ (the spectral counts distribution is that observed 
by construction). The required cold gas mass would be Mhi ~ 2 X 10^ -^©; from 
observations (Knapp et al. 1985) an upper limit of 1.5 X 10^ Mq for the HI mass has 
been estimated. Based on these results, we consider the hypothesis of an intrinsic 
absorption unlikely. 

6.2 X-ray Emission Sources in NGC4365 

A hot ISM doesn't need to be distributed as the stars; indeed, its Sx(-R) can 
be very flat, during the outflow phase (Ciotti et al. 1992). So, introducing the 
presence of hot gas (as discussed earlier, this gas would be hotter the the observed 
VCS, and therefore its presence is not required by the X-ray spectral analysis) 
could help reproduce the particular shape of the observed Sx(-R)- The amount and 
distribution of dark matter then play an important role, because - as simulations 
show - the gas density profile is flatter for larger and more extended (large /3 values) 
dark matter halos. 

The gas flow that allows the best agreement with both the observed X-ray 
spectrum and Sx(-R) must have the following characteristics: its Sx(-R) is mostly 
flat, and its Lx is not larger than ~ 2 X lO"^" erg s~-^, since a higher value cannot 
be reconciled with the observed spectrum. To reproduce the latter at high (> 1 
keV) and low (< 0.5 keV) energies, two further components must be added to the 
spectrum of the gas flow: a hard component, and a very soft component. The 
temperature and normalization of these additional components are determined by 
a best fit to the observed spectrum. 

After an accurate exploration of the (i2, /3, i^sn) parameter space, we found 
that the galaxy model best reproducing the data has [R = 9, /3 = 9.4, i^sN = 0.59). 
The X-ray properties of the three components are summarized in Table 4. Lx,vsc 
strongly depends on the assumed Nh value; here the value in Table 1 is assumed. 
The very soft emission can also be produced by a black body component of kTbi, = 50 
eV, rather than by a Raymond thermal component, and so the VSC temperature in 
this model is close to that of the supersoft sources (§5.1). To reproduce the whole 
Lx,vsc in this way we need ~ 200 supersoft sources. 

The resulting "best fit" spectrum is plotted in Fig. 16; x^/^ = 1-6, where 
= 18 is the number of degrees of freedom (the lower and higher energy spectral 
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bins have been excluded from the fit, see e.g. Fabbiano et al. 1994). The Sx(-R) 
resulting from the superposition of the three components is plotted in Fig. 17 
(x^/^ = 1-3, for V = 11); the spatial distribution of the 0.1 keV and the 5 keV 
components is the same as that of the stars. 

In conclusion, we suggest for this galaxy a SNIa's explosion rate of 0.6 the 
Tammann 1982 rate, and a mass model in which there is a considerable amount 
of dark matter, distributed mostly at large radii. For the reasons explained at 
§6.1, we stress that this result apply insofar the data support it, and needs more 
observational confirmation, such as can be obtained from the ROSAT HRI. The 
stellar velocity dispersion profile of NGC4365 is observed only within the central 
~ 30" (Bender and Surma 1992), so a useful check for the dark matter effects 
expected from our galaxy model is not possible. Bender and Surma 1992 also 
find that NGC4365 might have suffered a merging process in the past, based on 
the presence of a kinematically decoupled core component. If this is true, and a 
starburst followed the merger, the stellar population of NGC4365 might be peculiar. 

7. Conclusions 

In this paper we have investigated possible origins for the very soft X-ray 
emission [kT ~ 0.2 keV) recently discovered in low Lx/Lb early-type galaxies. We 
considered both an origin in gas flows, performing numerical simulations of their 
behavior tailored to two galaxies (NGC4697 and NGC4365), and an origin in stellar 
sources. Our findings are as follows: 

1) The X-ray luminosity of gas flows can be comparable to that of the very 
soft emission during the outflow or the partial wind dynamical phases. Both phases 
can be found at the present epoch only if the SNIa rate of explosion is lower than 
the standard (Tammann 1982) rate. The PW phase occurs only for values of the 
dark/luminous matter ratio up to 5-6, and the outflow for larger values. 

The average emission temperature of gas flows is comparable to that of the 
very soft emission only in galaxies with quite shallow potential wells, i.e. with low 
values of the central stellar velocity dispersion (cr^<200 km s~^). Due to the larger 
value, the average emission temperature of gas flows in NGC4365 (kT > 0.4 
keV) is higher than that of the very soft emission [kT < 0.3 keV at 99% confidence 
level) determined from the ROSAT pointed observation of NGC4365. So, by itself 
a hot ISM cannot generally explain the very soft emission. 

2) Concerning stellar sources, we investigated the possible contributions to the 
very soft emission coming from supersoft sources, late type stellar coronae, and RS 
CVn systems. 

The typical emission temperature of supersoft sources {kTbi, < 70 eV) is too 
low to explain entirely the very soft emission of early-type galaxies, although these 
objects might make an important contribution. In principle, if their nature is that 
of accreting white dwarfs, they could also account for the observed SNIa explosion 
rate. 
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M dwarf stars' coronae possess a thermal component at a temperature close 
to that of the very soft emission of early-type galaxies, but they are not likely to 
account entirely for it, because an age effect which reduces their X-ray luminosities 
is expected in old stellar populations. The initial mass function of M dwarfs should 
be much steeper than the Salpeter IMF to compensate for the age effect, but this 
may conflict with the observed M^/Lb ratios of ellipticals. 

RS CVn systems can also give a contribution, being the most X-ray luminous 
late type stellar sources, and possessing a soft [kT 0.2 keV) thermal component. 

3) We construct a detailed model for the X-ray emission of NGC4365, which 
can account for the results of the ROSAT PSPC pointed observation. We suggest 
the presence of three contributors, approximately equally important: evolved stellar 
sources such as low-mass X-ray binaries, producing the hard emission (energies > 1 
keV); a hot ISM, to better reproduce the X-ray surface brightness profile [kTgas ~ 
0.6 keV); stellar sources with very soft emission [kT 0.1 keV). A major fraction of 
the latter could come from supersoft sources such as those in 2). The galaxy model 
that produces a gas flow with the required properties has a dark-to-luminous mass 
ratio of 9, with the dark matter very broadly distributed, and a SNIa explosion rate 
of ~ 0.6 the Tammann rate. This galaxy model can be confirmed with ROSAT HRI 
pointed observations, which will give accurate measurements of the central Sx- 

4) In NGC4697, due to the low value of the central stellar velocity dispersion 
(cr^ = 186 km s~-^), the ISM is in a wind regime in most of the parameter space 
describing the dark matter presence and the supernova explosion rate. If we don't 
assume that considerable variations from galaxy to galaxy exist in the amount of 
dark matter and SNIa rate, the very soft emission should come entirely from stellar 
sources in NGC4697. This can be tested by future ROSAT pointed observations 
of this galaxy. 
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FIGURE CAPTIONS 



Figure 1: The Lx-Lb diagram for early-type galaxies; data (detections only) are 
from Fabbiano, Kim, and Trinchieri (1992). The dashed line is the luminosity 
expected from SNIa heating Lsn, assuming the Tammann (1982) rate of explosion. 
Group 1 galaxies are indicated with different symbols (full squares). 
Figure 2: The average spectral energy distribution of Group 1 galaxies, obtained 
by the Einstein IPC (open squares with errorbars). Also plotted are the two 
thermal components which best fit this distribution (dotted and dashed lines), and 
their superposition (solid line). 

Figure 3: 68% and 90% probability confidence contours for the temperatures (in 
keV) of the two components which best fit the average spectrum of Group 1 galaxies: 
the VSC (abscissa) and the hard component (ordinate). The increments above the 
Xmin value correspond to 3 interesting parameters, as explained in the text (Avni 
1976). 

Figure 4: a) as in Fig. 3, for the ROSAT PSPC spectral data of NGC4365; also 
the 99% probability confidence contour is shown, b) the ROSAT PSPC X-ray 
spectrum of NGC4365, together with the VSC (dotted line), the hard component 
(dashed line), and their superposition (solid line), at the best fit (from Fabbiano 
et al. 1994). 

Figure 5: The radial trend of the dark to luminous mass ratio Mh{< r)/M^[< r), 
for the cases R = 1 (lower lines) and R = 10 (upper lines); j3 = 2 for the dashed 
lines, and /3 = 10 for the dot-dashed lines. 

Figure 6: The (i2, i^sn) plane for NGC4697. The dashed lines enclose the region of 
PWs. The solid lines enclose the region where the models have Lx,gas ~ Lx,vsc] 
they are coincident in the outflow regime, when a small variation of i^sN causes a 
large variation of Lx, gas- 
Figure 7: The time evolution of the fundamental hydrodynamical quantities de- 
scribing the gas properties: velocity (v<0 for inflows, v>0 for outflows), temper- 
ature T and number density n. Plotted is the PW case [R = 5, i^sN = 0.27) for 
NGC4697. Dashed lines refer to i=0.75 Gyrs, dotted lines to i=4.5 Gyrs, and solid 
lines to the present epoch (i=15 Gyrs). 

Figure 8: The X-ray spectrum of the VSC from the analysis of Einstein spectral 
data (solid line), and that of a thermal component with kT = 0.5 keV, the 90% 
confidence upper limit on the VSC temperature (dotted line). The other lines are 
the gas spectra (convolved with the IPC response) for the following models for 
NGC4697: R = 9, i^sn = 0.38 (dashed line); R = Q, i?sN=0.28 (long dashed line); 
R = A, i^sN = 0.23 (dot-dashed line); R = 3, 'i?sN=0.2 (dot-long dashed line). In 
order to compare the shapes, each spectrum has been normalized to the number of 
counts in the channel where it has the maximum. 

Figure 9: The observed spectral count distribution for Group 1 galaxies (open 
squares with errorbars) compared with the superposition (solid line) of the con- 
volved spectral distributions of the gas flow in NGC4697, for a model with R = 3, 
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j3 = 2 and i^sN = 0.2 (dotted line), and of a hard [kT ~ 5 keV) thermal component 
(dashed line). 

Figure 10: The (i2, i^sn) plane for NGC4365. The solid lines enclose the region 
where the models have Lx,gas ~ Lx,vsc, and the dashed lines the PW region. 
Figure 11: The time evolution of the X-ray luminosity in the ROSAT sensitivity 
band for the following models for NGC4365: R=l, i^sN = 0.25 (sohd line); R=2, 
T?SN = 0.3 (dotted line); R = 3.5, i?sN = 0.36 (dashed line); R=9, i?sN = 0.68 (dot- 
dashed line). The three PWs can be recognized from the early cooling catastrophes 
(peaks in Lx)- To show the strong sensitivity to i^sN during the outflow phase, 
the case oi R = 5 has also been plotted (long dashed lines): in order of increasing 
Lx,gas-, the models have i^sN = 0.500, 0.498, 0.495 (this last is an inflow at the 
present epoch). The parallel solid lines represent the range oi Lx,vsc values allowed 
by the data. 

Figure 12: The observed spectral counts distribution for NGC4365 (open squares 
with errorbars) together with the superposition (solid line) of the distribution of a 
gas flow of the (i2 = 3, i^sN = 0.35) model (dotted line), and that of a hard compo- 
nent with kT = 5 keV (dashed line), both convolved with the PSPC instrumental 
response. In a) Nh = 1-6 X lO^" cm'^, and in b) Nh = 3.2 x 10^^ cm'^. 
Figure 13: 68%, 90%, and 99% probability confidence contours for the temper- 
atures of a hard thermal component (ordinate), and of a black body component 
(abscissa), derived fitting the ROSAT PSPC data of NGC4365 (see §5.1.1). 
Figure 14: The observed X-ray surface brightness profile of NGC4365 (open 
squares with errorbars), from Fabbiano et al. 1994; unbinned PI channels 26-247 
have been considered, to eliminate ghost effects. Also plotted are the surface bright- 
ness profiles of the hard (long dashed line) and soft (short dashed line) thermal 
components, derived from the best fit of the observed X-ray spectrum (§2), and 
their superposition (solid line). The thermal components are distributed as the 
optical profile of the galaxy (King 1978), and have been convolved with the PSPC 
response, appropriate for their respective spectral energy distribution. 
Figure 15: The same as in figure 14, but an intrinsic absorption corresponding to 
Nh = 5 X 10^^ cm~^ has been applied within the central 30" (see §6.1). 
Figure 16: The observed spectral energy distribution of NGC4365 (open squares 
with errorbars), together with the three components of the "best fit" model (§6.2) 
after convolution with the PSPC response: the hot ISM (dotted line), the very soft 
component [kT 0.1 keV, dot-dashed line), and the hard component [kT = 5 keV, 
dashed line); the solid line is their superposition. 

Figure 17: The X-ray surface brightness profile (solid line) of the superposition 
of the three components which best fit the observed X-ray properties of NGC4365 
(§6.2). Open squares are the observed counts, the long dashed line is the hard 
thermal component, the short dashed line is the very soft kT 0.1 keV thermal 
component, and the dotted line is the hot gas. 
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TABLE 1 



Observational Properties 



NAME 




„■ b 

) (km s" 


d Njf 
1) (Mpc) (1020cm-2) 


log 

(erg s"^) 


Lx.VSC f irp f 

(keV) 


NGC4697 
NGC4365 


1.95 
0.65 


186 
262 


37.4 2.1 
27 1.6 


40.97 
40.74 


0.15-0.50 <0.5 
0.3 - 0.7 0.14-0.21 



° total blue luminosity, from Fabbiano, Kim, and Trinchieri 1992 
central stellar velocity dispersion, from Whitmore, McElroy, and Tonry 1985 
distance, from Fabbiano, Kim, and Trinchieri 1992 

Galactic hydrogen column density in the line of sight of the galaxies, from 
Stark et al. 1992 

^ NGC4697 X-ray luminosity is the estimate derived from Einstein data by 
Fabbiano, Kim, and Trinchieri 1992, in the (0.2-4) keV range; Lx for NGC4365 
is derived from ROSAT data by Fabbiano et al. 1994, in the (0.1-2.0) keV range. 
From Einstein data, Fabbiano, Kim, and Trinchieri 1992 derived logix =40.42 for 
NGC4365, in the (0.2-4) keV range, assuming 1 keV thermal emission. 

^ 90% confidence range (see §2); the Galactic absorption in the line of sight of 
the galaxies has been assumed. 
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TABLE 2 



x(15 Gyrs), with i^sN = 1 







= 186 km 






= 262 km s 


-1 


R 





5 


10 





5 


10 


/5 














1 


0.12 


0.36 


0.60 


0.25 


0.71 


1.18 


6 


0.12 


0.31 


0.49 


0.25 


0.61 


0.98 


11 


0.12 


0.29 


0.46 


0.25 


0.58 


0.92 
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TABLE 3 



NGC4365: gas flows average emission temperat 
(t=15 Gyrs) 



R 




flow " 


hT 

"'-'-gas 






phase 


(keV) 


9 


0.68 





0.7 


5 


0.498 





0.5 


3 


0.35 


PW 


0.4 





0.18 


PW 


0.4 



outflow; PW= partial wind 
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TABLE 4 



X-ray emission components for NGC4365 





kT 
(keV) 


Lx{0.1 - 2 keV) 
(lO^" erg s-i) 


Lx/Lx,tot 


soft sources 


^ 0.1 


2.0 


0.35 


hard sources 


5 


2.1 


0.37 


hot ISM 


0.6 


1.6 


0.28 



Note: the X-ray luminosities, and the average X-ray emission temperature of 
the hot gas, are calculated considering the same region used by Fabbiano et al. 1993 
(a circle of 250" radius). 
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